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MILAT I’KODUCTIVl iY KSTIMiVFKS US I NO I^NNDSAT DATA 
TYl’i: I I I'RlHiUKSS UKI'OUT 
16 May 1977 - 15 (K tobt-r 1977 

Tlu* followiii}' ropori scivck as tho ninth Type 11 l*iop,rt'8s Ufpurt for 
l.aiulsat Follow-on 1 nvost 1 p.at ion /^2U62b whlili Is t'ntltlod "Wlioat 
t’rotluft I V 1 1 y Kslimatos llslnj; l.andsal Data". 

Tills Invost Igat Ion has sovi?riil ohjoct I vos , including tho fnllowini;: 

1. To dovelop t I'ohn l(pic*s and procfdu"t*s for using Landsat data to 
fstlinato cliaract i-r 1st i<‘s of wlioat oanopics which are correlated 
wltli potential wheat grain yield. 

2. To denionsl r.it e the usefulness of Landsat data for estimation of 
wlntet wheat yield: 

a. for Irrigated and for non-lrr 1 gat ed test sites 

b. for two different years with varying weather conditions. 

1.0 I’ROm.KMS 

No significant problems were encountered during this reporting period. 

2.0 SIGNIFICANT RKSULTS 

Large area Landsat yield estimates have been gene'.ated. These results 
have been compared with estimates computed using tlie NOAA Center for 
(H iinat ic and Knvlronmental Assessmcmt (CCEA) meteorological yield model, 
both of these estimates have heen compared witli Kansas Crop and Livestock 
Reporting Service (KCLKS) estimates of yield. In an attempt to assess 
the relative and absolute accuracy of the Landsat and CCF.A est Im.atc's. 

The results to date are inconclusive. 

A large area direct wheat prediction procedure has been implemented. 
Tills procedure potentially overcomes many of the serious problems (e.g., 
small fields and cloud cover over specific sites) being f.iced by other 
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ava i labli- ap|>ro.irlu*s. Iiiiti.il roKiillH have pro'Viced u wheat production 
estimate (dtnparahle witli the Kt'l.KS eHtImato. 

1.0 A(:Tivrnr:s 

Th»' aetlvitieH we were Involved In during this reporting period 
concerned two general topicn: 1) large area yield eHtlmaten; and 

2) rlirect haodsat wlieat production eHtiinati'B. The following material In 
a cursiiry description of tlioae activities. More details will he avallah''.* 
in t he I Inn 1 report . 

4.0 I.AHCi; ARICA YIKU) E.STIMVI’KS 

in order to make a comparison between l.andsat ef<t imati's of yi<*Id and 
met eoro !»>{’, lea 1 (Cff^Y) yield model est lmati>s of yield, we decldi'd to make 
botlt kinds of estimates t>ver a large area ftrr which a<ieqiiate "true" 
values of .average yield wirre available. We ciurse to use tlie Central Crop 
Iteporting District of Kansas as our basis for comparison. Yield <*st Imates 
for all counties In the (lentral CRD are available from the Kansas (.’rop and 
I.ivestock Reporting Service, and they were used as the "correct" v.iiues. 

In the following material, we will first discuss l.andsat estimates of 
yield over the Central CRD, and then CCEA meteorological yield model 
estimates of yield. The two methods will then be compared. 

4.1 1.ANDSAT LARGE AREA YIELD ESTIMATES 

The Central Crop Reporting District (CRD) of Kansas was chosen for a 
l.andsat large area yield prediction demonstration. One rea.son this area 
was chosen Is that, of the Kansas CRD's, It best satisfied the requirement 
for adequate "training" data. Information t)ii Individual field yle’d which 
is necess.iry In order to calibrate a l.andsat wheat yield relation was 
.'ivai l.'ible for three sltt^s within the Central CRD. 

It was decldt*d to carry out tin* yield prediction test using early M.iy 
1976 Landsat data. l.andsat yield pri'dlctlon on the test sites was based 
on a regre.sslon relation betVt?oen the l.andsat green measure, SQ7S, and 
farmers' combine weight estimates of yield per harvested acre on the 
training sites 
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Ki'cauHi* of tlx* |iuhkII)Ii> variation in oxtoinnl offocls hikIi oh 
atmospluT Ic Iwv.i avor tlio trainint’ and toKt nitos, it in |>oi>Hil>1c tlint a 
c’oi 1 1 ‘i't ion for Kucli factora woold bo rocjuirod. Accordingly, tlio procodnro 
of training and tiating a I-andHat yield alBorltiim over tlio Central CRi) waa 
lepeated it', ing data t liat waa corrected for amount of iiaze in tlie atmoapliere 
l»y a reci'iitly developi'd KRIM liaze normalizing program calied XSTAR. 

Wlien tile training algoritlinia were applied to tlie teat aitea, the 
yii'ld predict iona t’lat reaiilted are aiiown in Table 1. The l.andant 
eat imatea appear to be aenaitive to yield variation, alnce the unror- 
rected l.anda.it county average eatlraatea liave a variance of A.A9 compared 
witli a variance of 6.35 for KCI.RS eat imatea. Tlie Individual correlation 
lietween l.andaat yield eat imatea of a particular teat alte and county KCI.RS 
average yielda la not large for elflier the uncorrected (r ■ 0.25) l.andaat 
data or for tlie XSTAR corrected (r ■ 0.08) l.andaat data. 

It la not eaaential that these county eatlmates be hlglily correlated 
for tlie teclinique to be working, since a small sample in a county may not 
bi* repreaentat Ive of tiie wliole county. What la hoped, liowever, la tliat 
tlieae county samples, wlien appropriately aggregated, will be good 
Indicators of average yield over the entire Central CRD. In order to 
Invest igatc tills poaslblllty, tlie individual county yield eatlmates were 
weighted by tlie number of harvested acres of wheat for the respective 
county, and aggregated to determine an average value of yield for tlie 
Central CRD. 

The l.andaat average value of weighted county yields was then compared 
with tlie KCI.RS average yield, uaing a t-test . The hypothesia waa that the 
meana are identical. This hypothesis was barely accepted at the 5 percent 
level for the uncorrected l.andsat estimate of average yield. There appears 
to be a bi.ia in the l.andsat estimates of yield since most l.andsat estimates 
were too high uaing both XSTAR data (+A.2 bu/acre) and using uncorrected 
dat.i (+2.9 bu/acre). Apparently the source of Mas was not one that could 
be corrected by only accounting for atmospheric effecta (haze). 
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TABI.K 1 

KCLKS Actual YIcIcJh ami l.atulKat Predicted Yields (unweighted) 

Actual Yield Land sat Yield 


County 

_ (kci,rs)_ 

Uncorrected 

X_S^r/^ _Co r r cc t ed 

Sal Ine 

27.5 

37.8 

38.0 

nils 

30.6 

33.7 

36.6 

Marion 

29.3 

30.7 

31. 3 

Mcl’lierson 

28.5 

16.1 

36.3 

Kush 

30.8 

35.1 

32.5 

Rice 

36.3 

35.2 

37.8 

Russel 1 

36. 5 

32.9 

33.5 

K1 Isworth 

30.5 

32.9 

36.1 

(bounty Average 

30.8 

33.7 

36.7 

Standard Deviation 

2.5 

2.1 

2.5 
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4.2 ACIROMKThOKOl.OGlCAI. MOUi:i. Y I KU) ESTtM/MKS 

SIiu’o a^’ronift I'liroloj'lcn 1 ylt-lcl mcul«>lH arc bu f ri'qiit»nt ly used, llit>y 
ar«* ill Horao HfiiMf a yarclHLlck with wlilrh to ovaliiato al t iTiiat i ve a|>pruacii(*R. 
Ill t lu* fallowing si'ctloiiH wc* will di-scrliu’ tiic rosulls of impli'iiuMitat Ion 
of nil agromot ylold model, and wc wilt subsequently compare those results 
witli I.aiidsat results. 

The agroroeteorologlral yield model wlilrli was implemented for this 
projeet was a model developed to operate In Kansas by tlie Center for 
(!limatie. and Environmental Assesameiit (CCEA) of the National Oeeanogri'pble 
and Atniospberlc Administration*. The model was Implemented for the 
Central Crop Reporting District (CRD) of Kansas using all readily available 
data from meteorological stations scattered tlirougb the central CRD. We 
chose to implement the model for May truncation, since wc Intended to 
examine late April and early May Eandsat data, and since for April 
truncation no CCEA model was available. 

After the CCEA estimates were calculated for each meteorological 
station, an average value was obtained for each county with more than one 
meteorological station. 

It was possible to get coni])lete weather data from ten of the meteoro- 
logical stations located In the central CRD. CCEA agromet model estimates 
of yield were calculated for the ten meteorological stations and iv* .. ed 
with KCI.RS county estimates. Tlie unweighted CCE<\ estimates and the KCLKS 
estimates were found to have a non-s Ign I f leant correlation (r ■ 0.09). 

I.esn than 1 percent of the variance In KCLKS estimates Is accounted for 
by the agromet (CCl-^) estimates. 

The CCEA estimates are very stable, or conservative. The variance 
in county CCEA estimates is 1.01, wher»*as the variance In KCLKS county 
estimates Is 6.35. One might have exjiectcd the point samples (CCEA 
estimates geiu*rally from single meteorological stations) to be more 
variable than large area averages (KCLKS county estimates). However, the 


* NASA/JSC, 1975, "Wheat Yield Models for the United States", 
.1 .SC- 11 656. 


I.AClE-00431 , 
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CCKA <n'ronu‘t juTturb.it Ion niodol Ih not vory HciiHittvo to c'li;iii}*rH in 
wiMlhor. An nctd i t i onn 1 oxamplo of tliln rolut Ivo InHcnw 1 1 f v 1 1 y In that 
If thoro li.id hoi-n no procf pitat ion bc-twoon Auppint and Kidiruary, th** CCKA 
modol would havi’ prodlctod a yield reduction from normal yield of only 
i./ ho/aere. In reality, hucIi h situation would likely have had catn- 
Hlrophlc effects In yield. 

'’he individual county samp»e estimates of yield were subsequently 
weighted by tin wheat ncreagt* harvested (In 1976) in the county corre- 
sponding with the meteorological station(s). The estimates were then 
aggregated to .i single estimate for the central CKI), as was done using 
l.andsat estimates. Despite tin* apparent insensitivity of the CCKA model 
to meteorological variations (or perhaps because of It), and despite the 
low correlation bt'tween CCKA <ind KCKKS estimates, tin* aver ige weighted 
CCKA value of yield is not far removed from the KCLRS estimate. The 
difference is 1.6 bu/acre, which has a I’-value of 0.18. Therefore, we 
accept the estimate of yield as being not statistically significantly 
different from yield. 

The above discussion indicates both tin? advantages and d 1 sidvantagi*s 
of an agromet perturbation model of the type implemented. Its stability 
and relative freedom from a constant bias generally guarantees that it 
will not be far in error in reas(*nahly "nornal" years. However, its 
conservativeness also precludes it from adetpiately reflecting the effects 
on yield of large deviations from normal weather. 

4.3 COMPARISON OF ACROMKT AND LANDSAT KSTl MATES 

The preceding discussion (Sections 4.1, 4.2) may not furnish us with 
definitive answers that reflect the general relative utility of agromet 
and l.andsat yield estimates. For example, whether individual county 
estimated of yield using the two techniques are correlated with KCLRS 
yield or not may not he terribly relevant, because of the sampling 
schemes used. Similarly, the accuracy of prediction of weighted average 
yield is not necessarily definitive. This Is due to the fact that moat 


6 



NkT • I A Ml.N •«« tC lk<t » TMf ti»Mv|li%tt f > < M< 


'^rRiM 


of itio " Informal Ion" In tills partlciilnr t»>Ht s<u-ms to be In tlie acrcaKf 
wolf'll! ing fai'loi H, wliirli li.ivo substanl Ini ly lary,or coof f Ic lont of 
V. trial ion than do ilu* ylohl oHlIinatcs (Tablo 2). TluTcforo, tlu* county 
with tlu* larj'oHt liarveslod wheat acreai'c tends to have the lar^'est 
weltihted yield esllmale, rey.ardleHs ol the type of yield (nitlmate (handsnt, 
KCI.RS, or CCKA), and conversely for the county with the smallest harvested 
wheat acreage. This situation results In, for example, unweighted CCKA 
estimates luivlng a correlation with KCl.KI) estimates of 0.09, and the 
corresponding wi'lghted estimates having a correlation of 0.92. 

Despite these difficulties In Interpretation, the results do shed 
some light on characlt r 1st Ics of the two approaches that might be fairly 
general In nature. Spev i f j ca 1 I y , the agromet model is characterized by 
relative lack of consistent yield t*rror (bias) and by insei sit Ivlt y to 
large changes In yield. The present lainds.it model Is characterized by 
potentially large yield bias and by high sensitivity to changes In yield. 

In other words, either approach has advantages and disadvantages. Klther 
approach might be modified to reduce Its disadvantages. Also, agromet and 
I.andsat Information could be used together to estimate yield. These 
possibl lit lea are being further explored, and will be documented In the 
final report. 

S.O DIRECT LANDSAT WHEAT I’RDDUCTION ESTIMATES 

Thus far we have discussed only the ability to forecast wheat yield 
(per acre) using I.andsat data. By Itself, this Information would be 
valuable as part of a system for forecasting wheat production. However, 
our work to this point has suggested a method for utilizing the relation- 
ship between I.andsat data ami yield, together with other relationships, 
to effect direct I.andsat forecasts of winter wheat production which may 
overcome certain troublesome problems In some of the existing approaches. 

The existing approaches tend to separate the task of forecasting 
Into two separate subsystems consisting of; (1) wheat acreage determi- 
nation; and (2) regional average determination of per acre yield. The 
ajiproach discussed below could make It possible to determine production 
on a plxel-by-p?.xel basis, using enrly-season I.andsat data, with a single 
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TABLK 2 

Coefficient of Variation (o/m) for Product Ion-Related ParumeterH 
from Count li‘H Wltliln tlie Central CRD 


Parameter 

o/m 

KCI.RS 

0.08 

CCKA 

0.03 

LANDSAT 

0.06 

Aereaj'e 

0.22 
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procfriHlng step. TIiuh, It may l>»*corai* poHaiblt* to Hiirvi-y JarKc urenH 
Hiirli aa n atati' or country much more ' M>nomlca1 ly than at prcHcnt, and 
acliicvo more timely information. Wlwu followa ia a dtacuaslon of tltc 
rationale of t lie au)>p.eated approach, and a d(>monatrat ion of ita Initial 
implcraontat ion. 

One of tiu* ideaa heliind tlie direct wtieat product ion approach using 
l.anduat data la tliat an appropriate value of production can l»c determined 
for eacli pixel in tlio act*ne, perliapii witliout even the m-ed to apeclfly 
wlietlier tile pixel la wlieat . 

We liave prevloualy aliown tliat several I.andaat transforms are good 
meaaurt'a of green vi'gelat Ive covc'r, and tliat cover In turn is strongly 
relaU'd to wlieat yield. Olven the knowledge of the area covered by a 
pixel the estimate of yield on a per pixel hasla can be directly converted 
to production. An additional f.ict la that in winter wheat regions such as 
Kanaaa, wheat tends to develop aigntf leant green cover sooner than most 
non-wheat fields and can therefore be eaally d 1 at i nguiahed . (Wlieat 
c laasif Icat Ion accuracies of 92 and percent were achieved on two 
Kansas sites using only the l.andsat SQ75 green measure.) Thus, If a 
production-predictive relation (developed on wheat fields) Is applit^d 
to non-wheat pixels, a very low production indication would be ^.xpected, 
and might be a negligible source of error. If applied to pixels falling 
cn ii bound.iry betw(*en whi*at and non-wheat, an appropriate intermi’d lat e 
value of green cover, and thus, intermediate average production would be 
estimated. This Intermi*diate value of production could approximate the 
total amount of wheat production represented by the pixel, which covers 
an are.a only p;irtliilly planted to whe.'it. Thus, pixeln would tend to 
contribute only their fair share of the total production estimate. 

Ah a part of this procedure tt is necessary to establish the production 
predictive relationship on an area wliere ground truth information is 
available*. With the relationship established, the present approach is to 


* In an operational env Ironnuuit , several carefully selected sites and data 
from previous years should satisfy the need for training. 
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Nflt‘ct a thr(>Hlu*ld below wtiich no wlieal iirculiK'l ion in JiHniKneil to a t’.^ven 
pixel. The need for nuch n tltreiihold 1 h dictated by tlie fjict that, in 
p.enernl, Home non-wlieat pixelH Renerafe l.andHat green me^ aren wbicti fall 
above tlu»HC of Home lt*w production wiu-at pixelH. The tlirc’Hbold value 1 h 
H eiecled to r<iuBe err<»rH of omiHKion and core l•.8lon to c«>mpenHate. 

Ah an Initial teat of tlie direct production forerant procedure, tlie 
above npproacli wan employed uHing tl>e SQ75 green measure t»n a portion of 
the 6 May 1976 l.andHat data for Slt«* A. Kmploylng the reKultIng relot lon- 
hbip on all of Site A a production forecast of 42,700 busbelH was m.ide. 

Tlila compared favorably wltli the actual production of 40,600 bushelH for 
thiH nltc, an error of only h.2 prr'ent. In addition we applied the name 
procedure to a different Hite (fUte R) ag.iln uHing 6 May 197b l.andHat data. 
The reHultlng production eHtlm.iteu for these testn are Hliovm in Table* 1. 
Note that the total production estlmatc'd for the two te«tH with separate 
training was within 1.6 percent of the correct total product l«»n, well with- 
in I.ACIK desired accuracy*. 

A further teat of the handnat direct wheat product ion approach was 
performed over ten of the eleven counties of the entire central CRD** 

UHing mid-April data. In this case, training was accomplished using six 
Hites, and the "toHt" data Het conslntvd of a 7.7 percent sample of the 
avall.'ible l.endsat pixels over the ten countie « of the central CUD***. 

The resulting l.andHat production estimate for the ten counties was within 
2.6 percent of the final KCI.RS estimate. TIiIh error is leaf than the 
Hpread of the several preliminary KCI.RS estimates made In the months 
fol lowing Apr 1 1 . 


* MacDonald, R. B. , Hall, K. , and Krb, R. B. , "The Large Area Crop 

Inventory Kxper iment (I.ACIK) - An Assessment After One Year of Operation", 
l’roce»-d ings of Tenth International Symposium on Remote Sensing of 
Knvlronment, Environmental Research InHtltute of Michigan, Ann Arbor, 
Michigan, 197^. 

** No Landsat data was available for the Sout heast ern -m<»st county. 

*** Further detail •« will be given in the final report. 
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TABLE 3 

Iiiilial RchuU from KRIM Direct Wlieat Production ForecaHl Procedure 
(Two LACIK IntenHlvo TohI Sltea) 



IJ^NDSAT 

OverpaHH 

True 

Product Ion 

ERIM 

Production 

Kort^aat^^ 

Error 

JAL- 

A 

6 May 76 

40,600 bu 

42,700 bu 

5.2 

B 

6 May 76 

27,900 bu 

24,700 bu 

11.5 

A+B 

6 May 76 

6H,:‘>00 bu 

67,400 bu 

1.6 
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The* ahnvi* ti'HlH of the* KKIM l.aiulHat whoat production <*st im.it Ion 
Kyatom arc certainly not definitive. Mary more tcHta in dlff**rent 
HltnatioiiK need to lx* carried out in order to aHHeitH tin' conHlatenry of 
performance. However, the preliminary indicat ions b.ased on onr limited 
te'rttH Ktvo encouragement tli.it tiie direct whe.it product ion npproacli 
uHing early-ceaHun l.andaat data la worth puraulng. 

(>.0 FirrUHK 1M.ANS 

Technical efforta of tiiia project are in the final Htugea. Tln^ae 
effortn will be fully documented in a final report, wliicli ia currently 
being prepared. 

7.0 PU Bl.lCAT I ON'S/PKKSKNIAT I ONS 

Mr. R. K. pka and Dr. J. K. Colwell attended tlie l.aiulHat Crop 

Condition aai Yield Briefing iield at NASA lleadquar tera on Septenilier 27, 
I'l-V. Mr. Nalepkn pn*Hentod a review of our actlvitlea entitled "World- 
wide Wlieat Production Korecaata Uaing L.nulaal Data". 
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